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Analytical calculation of electron’s group velocity surfaces in uniform strained graphene 
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Electron group velocity for graphene under uniform strain is obtained analitically by using the Tight-Binding 
approximation. Such closed analytical expressions are useful in order to calculate electronic, thermal and optical 
properties of strained graphene. These results allow to understand the behavior of electrons when graphene is 
subjected to strong strain and nonlinear corrections, for which the usual Dirac approach is not longer valid. 
Some particular cases of uniaxial and shear strain were analized. The evolution of the electron group velocity 
indicates a break up of the trigonal warping symmetry, which is replaced by a warping consistent with the 
symmetry of the strained reciprocal lattice. The Fermi velocity becomes strongly anisotropic, i.e, for a strong 
pure shear-strain (20% of the lattice parameter), the two inequivalent Dirac cones merge and the Fermi velocity 
is zero in one of the principal axis of deformation. We found that non-linear terms are essential to describe the 
effects of deformation for electrons near or at the Fermi energy. 

PACS numbers: 


INTRODUCTION 

Graphene was the first two-dimensional (2D) crystal dis¬ 
covered m. It has been broadly studied due to the observed 
peculiar physical properties GHS- The electronic proper¬ 
ties are mainly determined by electrons at the Fermi energy 
(4). For graphene, such electrons have momentum near or at 
the high symmetry points of Brillouin Zone (BZ). This be- 
havoir can be modeled by a Dirac Hamiltonian 0, where 
electrons behave as massless Dirac fermions with a Fermi ve¬ 
locity Vf ~ 1 x 10 6 m/s, which plays the role of the speed of 
light. In this model, the Fermi velocity is a constant parameter. 
However, this is not longer true when graphene has corruga¬ 
tions (curved space) or is stretched, since these deformations 
give rise to a space-dependent Fermi velocity 0 , suggesting 
changes in the electronic conductivity. Furthermore, in the 
case of stretching, a bandgap opening is observed 00 . Such 
results open the possibility for doing “strain engineering” in 
order to tailor the electronic properties and thus control the 
electron transport mm. 

Several theoretical approaches have been proposed to study 
deformations in graphene 16 HT6H2TI . The most common one 
is a combination of the Tight-Binding Hamiltonian (TB) and 
linear elasticity to derive a Dirac effective equation 0. Under 
such approach, pseudomagnetic fields appear, although lattice 
deformations were not included in the original derivation lt22l . 
In the case of strain, recent works have included these con¬ 
siderations starting from different treatments mmmm. 
Still, there are some problems with such approach l24l be¬ 
cause a common confusion is the assumption that the Dirac 
cone tips Kp> in the new deformed lattice coincides with the 
strained high symmetry points K and K'. 

In the present work, we calculate analitically the group ve¬ 
locity surfaces for the important case of uniform strain, which 
can be solved without the usual perturbative analysis of the 
Dirac equation. For this goal, the TB approximation has been 
used. The Fermi velocities are obtained by looking at the ap¬ 


preciate points in the reciprocal space. Thus, our results are 
more general and include the Dirac theory of strain as a limit¬ 
ing case. Additionally, we found that even for a realistic value 
of pure shear strain, a mixed Dirac-Schroedinger behavior can 
arise, suggesting that the Dirac theory has to be modified. In 
fact, this behavior has been obtained in other cases 01251. 

The layout of this paper is the following. In Sec. 2, we 
describe electron behavior in graphene under uniform strain. 
Then a dispersion relation is obtained using TB aproxima- 
tions. For this relation, we display the surfaces and contour 
plots for the particular case of pure shear strain. In Sec. 3, we 
derive the group velocity for electrons, and analyze the pure 
uniaxial and pure shear strain cases. Finally, in Sec. 4, we 
give our conclusions. 


ELECTRONS IN STRAINED GRAPHENE 

Graphene is formed by a single atomic layer of carbon 
atoms arranged in an hexagonal estructure, as shown in Fig. 
|T| The structure can be described in terms of two triangular 
sub-lattices, A and B, with a basis of two atoms per unit cell. 
The lattice unit vectors are given by ai and a 2 and the three 
nearest-neighbor vectors can be written as <5j, S 2 and S 3 0. 
Likewise, there are two reciprocal-lattice vectors given by 
and b>, and two inequivalent special points at the corners of 
the graphene BZ, called high symmetry points Kq and Kq 
0- For unstrained graphene, the tips of the Dirac cones (or 
the Dirac points Kp>) are located at the Kq and K' 0 points. 

In the case of a uniform strain, if the vector r represents the 
positions of the carbon atoms in the undeformed graphene, its 
deformed counterpart is given by r 1 = (I + e) ■ r, where I is 
the 2x2 identity matrix and e is the uniform strain tensor. The 
lattice unit and nearest-neighbor vectors are thus as follows, 
a[ = (I + e) • a t (i = 1, 2) and S' n = (I + e) • S n (n = 
1,2,3), while the reciprocal-lattice vectors are deformed as, 
b' = (I + e)” 1 ■ bi (i = 1,2). The new high symmetry 
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FIG. 1: (Color online) Graphene lattice and the sublattices A and B. 
The associated unitary and first neighbour vectors are also shown. 


points in the comers of the first BZ of the strained reciprocal 
lattice are obtained by construction of the Wigner-Seitz prim¬ 
itive cell and can written in general as 

K = M~ 1 C 1 and K' = M~ 1 C 2 (1) 


with 


and 


Mi = 


(&')* 

(6i)x + (^ 2 )® 


0 K)v 

(b'l)y + (b '2 )y 



( IMII 

m+by 


where (b[) x and (6')^ are the x and y components of the de¬ 
formed reciprocal vectors 6' (i = 1,2). 

To obtain the electronic properties of graphene under uni¬ 
form strain, we use the nearest-neighbor TB Hamiltonian l24l 


H = -Y,t' n albr' +S ' n + H.c, (2) 

r' ,n 


where r' runs over all sites of the deformed honeycomb lattice 
and the hopping integral t' n varies due to the modification in 
the carbon-carbon distances as t' n = t exp[— f}(\8' n \/a — 1)] 
li26l . with /3 « 3 and £ ss 2.7 eV is the hopping energy for 
unstrained graphene O. The operators n' r , and b r ' + g> corre¬ 
spond to creating and annihilating electrons on the A sublat¬ 
tice position r’ and B sublattice position r 1 + S' n , respectively. 
Now, using the Fourier representation for these operators, the 
previous Hamiltonian can be written as El 

H = -J2 t' n e- ik < i+ ^a{b k + H.c., (3) 

fe,n 


which finally leads to the closed dispersion relation for 
graphene under uniform strain El 


E(k) = ± 




-ifc-(J+e)-<5„ 


(4) 


It can be stressed that in our work, the latest equation has been 
developed more explicitly, which leads to 

E(k) = ±^ + g(k), (5) 


where 


3 3 

g( k ) = ^ cos [ fc • i 1 + e) • (8n - <5 S )] 

n—1s>n 


and 


7 = XX 2 . 

n—1 

Let us now explore the strain effects on the energy dispersion 
relation. As it was previously explained in the introduction, 
when an uniform strain is applied the reciprocal lattice is also 
strained. Thus, the first BZ is modified, i.e., its original hexag¬ 
onal form is varied to a polygonal form, as shown for the par¬ 
ticular cases of pure shear strain along the armchair direction: 

= £yy = 0, e xy = e yx = 0.1 (see red lines Fig. |2]a)) and 
£xy = £yx = 0.2 (see red lines Fig. [2]b)). In the same figure 
[2] along with the first BZ, we present the contour plot of the 
energy obtained from Eq. 0- 

Once the first BZ and the energy surfaces are obtained, 
we need to locate the Dirac points Kp using the condition 
E{Kd) = Ep, which corresponds to electrons at the Fermi 
level. By applying this condition, the Dirac points are indi¬ 
cated as pink circles in Fig. [2] The most important conclusion 
from the figure is that such points are no longer located at the 
high symmetry points K and K' (red circles) of the corners 
of the first BZ (Eq. |T|, since they are shifted to the saddle 
point. This plot ilustrates an issue that has been disregarded 
in the literature concerning graphene. 

As deformation increases (up to 20%), the Dirac points 
merge into the saddle point and a gap opens, which is con¬ 
sistent with the results in references 0E3- Furthermore, in 
this critical point, the dispersion relation is linear along one 
direction ((relativistic Dirac behavior) and quadratic along the 
other one (nonrelativistic Schrodinger behavior). Therefore, 
the Dirac theory needs to be modified. 

Summarizing the above, the effects caused in graphene un¬ 
der uniform strain are the following: 

• i) The Dirac points are shifted from the corners of the 
strained BZ. 

• ii) The Dirac equation is no longer suitable for long 
strain (> 20%), since for particular cases a Dirac- 
Schrodinger behavior is observed and furthermore one 
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FIG. 2: (Color online) Contour plot for the energy for a shear strain given by t xx = t yy = 0, and a) t xy =0.1 and b) t xy = 0.2. The first 
Brioullin zone of the strained reciprocal lattice is presented with red lines. The high symmetry points K and K' (Eq. [lj are indicated with red 
circles. The pink circles correspond to the position of the Dirac cones Kd, where the Fermi energy is located. Notice how two Dirac cones 
merges into one in case b) and do not have the same position as K and K'. 


might expect significant nonlinear corrections. It fol¬ 
lows that the anisotropic Fermi velocity is not longer 
valid in these regimens. Therefore, we must consider a 
more general velocity to understand the electron behav¬ 
ior. This is done through the calculation of the group 
velocity, as we will discuss in the following section. 

GROUP VELOCITY 

In the literature, the basic properties of electron transport 
phenomena in a crystal are described in terms of Bloch waves 
with wave vectors k l28l . Using these waves, we can build a 
dispersive wave packet with a certain group velocity. It can be 
shown quite generally that the mean electron velocity is given 
by the group velocity of the wave packet t28l 

v{k) = V k E{k), (6) 

where V& is the gradient operator in fc-space. From this equa¬ 
tion, the real-space motion of the electron can be described. 
Here we are interested in the behavior of electrons in graphene 
under uniform strain. Thus, by substituting Eq. ([5| into Eq. 
(|6) we obtain 

V(fc)=± 2 ^fe) Vfc5(fe) ' (?) 

The components x and y of v{k) are given by 

3 3 

vi(k) =±£Et(l + e u )(S l n - S l s ) + e lm (Sn ~ O 

n =1 s>n 

w J J Sin l fe • i 1 + e ) • (<5n - <*«)] 

X T n l s / -;- TT~\ 


where l,m = {a;, y} and l ^ to. The group velocity norm is 
given by v(k) = |u(fe)|. 

In the Fig. [3]we plot the surfaces and contour of the group 
velocity norm v{k) for pure graphene. It is important to note 
that at low energies and in the vicinity of the Dirac points, 
v(k) is isotropic, coinciding with the Fermi velocity. How¬ 
ever, as we move away from the Dirac point (corresponding 
to nonlinear corrections to the Dirac cone) a trigonal warping 
appears, giving rise to an anisotropic behavior (see contour 
Fig. 0- Furthermore, it can be observed that in the directions 
where the trigonal warping appear, v(k) increases, while in 
other directions decreases drastically. These results do not ap¬ 
pear when the Dirac theory is used. Therefore, if we want a 
more complete understanding of the behavior of electrons in 
the energy bands, nonlinear corrections and directions should 
be considered, since strain effects enhance such features, as 
we discuss below. 

We analyze the particular cases of pure uniaxial and pure 
shear strain with e xy = e yy = 0, e xx = 0.1,0.2 and 
e xx = £yy = 0, e xy = 0.1,0.2, as shown in Fig. [4] and 
Fig. 0 respectively. In these Fig. [4] and Fig. [5j the contours 
plots of the velocity norm v(k) are presented. Overimposed 
to these contours, we present the first BZ of the strained re¬ 
ciprocal lattice with red lines. Likewise, the high symmetry 
points (Eq. [TJ are indicated with red circles and the pink cir¬ 
cles correspond to the position of the Dirac point Kd where 
the Fermi energy is located. Notice again how the Kd points 
do not have the same position as the high symmetry points K 
and K'. The effect for shear strain is much more pronounced. 

On the other hand, the effects caused by the deformation in 
the velocity surfaces are the following: 

• i) As seen in [4] and Fig. [5] the Fermi velocity is no 
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FIG. 3: (Color online) Surface and contour plot of the group velocity norm for pure graphene. 




FIG. 4: (Color online) Countour plot of the group velocity for a uniaxial uniform strain given by e xy = e yy = 0, and with a) t xx = 0.1 and b) 
e xx = 0.2. The first Brioullin zone of the strained reciprocal lattice is presented with red lines. The high symmetry points K and K' (Eq. [lj 
are indicated with red circles. The pink circles correspond to the position of the Dirac cones Kd, where the Fermi energy is located. Notice 
how Kd do not have the same position as K and K'. 


longer isotropic. Instead, it becomes strongly anisotrop- 
ical. 

• ii) The surfaces do not display the trigonal symmetry, 
and instead, they present the symmetry of the corre¬ 
sponding strained reciprocal lattice. 

• iii) The trigonal warping observed in pure graphene dis¬ 
appears, as in (Fig. |4]a)). It reappears in (Fig. |4]b)) but 
with deformed angles which follow the symmetry of the 
strained reciprocal lattice. This new warping is strongly 
modulated, and it touches the Dirac points (Fig. [4] b)). 
Since the warping is associated with non-linear terms, 
this suggests that non-linearity is important in order to 


describe such cases. As a result, a pure Dirac equation 
is not longer valid. For the shear strain, in Fig. [5] a) the 
trigonal warping is deformed, until in Fig. [5]b) becomes 
two lobules reflecting the symmetry of the strained re¬ 
ciprocal space. 


• iv) When the Dirac cones merges by shear strain as in 
Fig. [5] b), the group velocity is zero along one of the 
pricipal axis of the deformation. This is just the conse¬ 
quence of the energy having a parabolic (Schrodinger) 
behavior with a gap opening (7J. 
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FIG. 5: (Color online) Countour plot of the group velocity for a shear strain given by t xx = 0, e yy = 0, and with a) t xy = 0.1 and b) 
t xy — 0.2. The first Brioullin zone of the strained reciprocal lattice is presented with red lines. The high symmetry points K and K' (Eq. [lj 
are indicated with red circles. The pink circles correspond to the position of the Dirac cones Kd, where the Fermi energy is located. Notice 
how two Dirac cones merge in case b) and do not have the same position as K and K'. From the countor plot, is clear that the Fermi velocity 
is constant (Dirac behavior) in the principal axis of the the shear, while it follows a parabolic behavior (Schroedinger) in the perpendicular 
direction. 


CONCLUSIONS 

We have obtained the electron’s group velocity for 
graphene under uniform strain using the tight-binding approx¬ 
imation. Our results indicate that the velocity is strongly 
anisotropic and that the trigonal warping is deformed to fol¬ 
low the symmetry of the strained reciprocal lattice. As strain 
increases, this warping touches the Dirac points. Thus, we 
found that non-linearity is very important in order to describe 
electrons in a proper way near the Fermi energy for strain, 
since the trigonal warping observed in graphene touches the 
Dirac point and gets modulated by the symmetry of the 
strained reciprocal lattice. As a result, a Dirac equation kind 
of approximation is not longer valid for such cases. 

Finally, our closed analytical expressions for the electron 
velocities are useful in order to calculate electronic, thermal 
and optical properties of strained graphene. 
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